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We study the effect of a constant uniform magnetic field on an electrically charged massive particle 
(an electron) bound by a potential well, which is described by means of a single attractive A5(r) 
potential. A transcendental equation that determines the electron energy spectrum is derived and 
solved. The electron wave function in the ground (bound) state is approximately constructed in a 
remarkable simple form. It is shown that there arises the probability current in the bound state in the 
presence of a uniform constant magnetic field. This (electric) current, being by the gauge invariant 
- - - quantity, must be observable and involve (and exercise infiuence on) the electron scattering. The 

probability current density resembles a stack of "pancake" vortices" whose circulating "currents" 
around the magnetic field direction (a-axes) are mostly confined within the plane 2 = 0. We also 
compute the tunnelling probability of electron from the bound to free state under a weak constant 
homogeneous electric field, which is parallel to the magnetic field. The model under consideration 
M , is briefiy discussed in two spatial dimensions. 
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I. INTRODUCTION 



The behavior of quantum nonrelativistic systems in external electromagnetic fields has attracted per- 
04 I manent interest in view of possible applications of the corresponding models in many phenomena of 

. quantum mechanics. Bound electron states play an important role in quantum systems in condensed 

matter. When the external field configuration has the cylindrical symmetry a natural assumption is that 
. the relevant quantum mechanical system is invariant along the symmetry (for example, z) axis and the 

' quantum mechanical problem then becomes essentially two-dimensional in the xy plane. Such is the 

, case when the external field configuration is a superposition of a constant uniform magnetic field and a 

' cylindrically symmetric potential. Nonrelativistic electrons in such external field backgrounds are good 

I quantum mechanical models for studying remarkable macroscopic quantum phenomena such as the frac- 

tional quantum Hall effect^ and high-temperature superconductivity^. A related problem is the behavior 
^ i' of electrons in trapping potentials in the presence of a constant uniform magnetic field in an effective 

mass approximation. Magnetic fields seem to be likely to effect on weakly bound electrons into singular 
' potentials of defects in the defect films'^ and solids'*. 

, Pure two-dimensional models also are of significant interest. The effect of magnetic fields on a weakly 

^ i' bound electron in two spatial dimensions was studied by us in Ref. 5. This model is of interest because it 

gives a good example of a nonrelativistic analog of the so-called dimensional transmutation phenomenon 
first discovered by S. Coleman and E. Weinberg in Ref. ^ in the massless scalar electrodynamics. Whereas 
above nonrelativistic systems can be described by the Schrodinger equation, relativistic systems, related 
to the Dirac Hamiltonian in 2-1-1 dimensions in a constant uniform magnetic field, show up in a certain 
' type of doped two-dimensional semimetals^. Similar problems are also related to a number of problems in 

quantum theory, for example, the parity violation, the theory of anyons (particles satisfying a fractional 
statistics), the Aharonov-Bohm effect^, and other. At last recently, a new type of spectral problem has 
been found^'*" in quantum mechanics of planar electrons in a superposition of constant uniform magnetic 
and cylindrically symmetric potential fields. For this new class of spectral problems, the so-called quasi- 
exactly solvable (QES) models the energy spectra exist if a certain relation between the parameters 
characterizing the intensity of the interaction of an electron with external fields hold and then solutions 
of the corresponding equations of quantum mechanics can be expressed as the product of a weight function 
and a finite polynomial. The some physical examples of QES models, which include the two-dimensional 
Schrodinger or Dirac equation for an electron in a superposition of a constant uniform magnetic field and 
an attractive Coulomb field, were studied in Ref. [ll| . 

The purpose of this paper is to study the effect of a constant uniform magnetic field on an electron 
bound by a single attractive X6{r) potential. We derive a simple transcendental equation determining the 
electron energy spectrum and construct the approximate wave function for a bound electron state in the 
presence of a constant uniform magnetic field. We show that the sizes of the electron localization region 
change and the probability current arises even when the electron is in the bound state in a superposition 
of a constant uniform magnetic field and a single attractive A5(r) potential. The probability current in 
three-dimensional space resembles a stack of "pancake vortices" whose circulating (around the z-axes) 
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"currents" are mostly confined in the weak magnetic field within the plane z — 0. 

The equation for determining the energy levels of the electron states is also obtained for the model 
under study in two spatial dimensions. We show that in difference from the three-dimensional case the 
binding energy is not analytical in A. 

The tunnelling of electron from the bound to free state under a weak electric field is of importance. 
In present paper we compute the tunnelling probability of electron from the bound to free state under a 
weak constant homogeneous electric field, which is parallel to the magnetic field. 



II. ELECTRON IN AN POTENTIAL WELL IN THE PRESENCE OF A CONSTANT 

UNIFORM MAGNETIC FIELD 

Let us consider an electron with the charge e < in an attractive singular potential of the form 

Uir) = -^^Sir) (1) 

and a constant uniform magnetic field B, which is specified in Cartesian coordinates as 

B = (0, 0, B) = V X A, A = {-yB, 0, 0) (2) 

In ([T]) A is a positive coupling constant of the length dimension, S{r) is the three-dimensional Dirac delta 
function, m > is the effective mass of an electron. It is well to note that the attractive XS{r) potential 
can be considered as the limit of sequence of appropriate narrow rectangular potential wells 

U{r) -Uo, r <R, U{r) =0, r > R, (3) 

then the parameter A is expressed via Uo and R as 

_ 2mUoR^ 
A — Q . 



The Schrodinger equation is 

The electron wave function in magnetic field ([2]) can be found in the form^^ 



1 

2m 



(5) 



where 



E^.^n.(n + \y§- + sn.£- (6) 

is the energy eigenvalue, co = \eB\/mc, me is the mass of free electron, pi,p3 are the eigenvalues of 
generalized momentum operator and s — ±1 is the spin quantum number. Note that for nif. — m and for 
signeB < 0, all the energy levels except one with n = 0, s = — 1 are doubly degenerate: the coincident 
levels are those with n s = I and n + 1, s = — 1. In given paper we consider the case mg = m. Then, 
the energy eigenvalues depend only on the number n. 

Note that pi is constrained by |pi| < eBL/c (see, Ref. [12]). The functions 

are expressed through the Hermite polynomials _ff„(z), the integer n — 0, 1, 2, . . . indicates the Landau 
level number, 



he h 
JeBl " V mw 
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is the so-called magnetic length and j/o = —cpi/eB. It should be reminded that the classical trajectory 
of electron in the xy plane is a circle. The quantity t/o corresponds to the classical y - coordinate of the 
circle center. All the electron states ([5]) are not localized in the x, z-directions. 
Solutions of Eq. (g]) are sought in the form 

oc „ 

*£;(r) = X! / dpidp3CEnp->JJnpir) = ^ C£;„pV„p(r), (7) 

n— n,p 

where ^„p(r) is the spatial part of wave functions ([5]). 
Then, for coefficients Ceup, one obtains 



CEnpin + cpl + b)^XoJ2 CEikViiO)VniO), (8) 

where 



Lk 



^ h = -^, Ao = , / , , vm = Vi{y = 0). 



Note that terms on the right of Eq. ([5]) are the matrix elements of the interaction operator of electron with 
the singular potential H]). Correctly these matrix elements must be obtained at first for the interaction 
potential ([3]) with finite Uq, R and then by proceeding to the limits R —^ 0, Uq oo. But for finite 
Uq and R the main contribution in the matrix element Unp^^n'p'^ gives the quantum number regions 
|p3 — Pgji? < \/n — n'R < a because the matrix element falls off quickly due to the oscillations outside 
the above regions. Moreover, for the sharp potentials, one can obtain Unp^^n'p'^ ~ ~A for n,n' < N and 
P3,P3 < P and Unp^^n'p'^ — for n,n' > N and P3,P3 > P, X ^ N, P oo as i? ^ 0. So, for such 
potentials the sum taken over n on the right of Eq. ^ is limited by the value N and the bound state 
energy becomes depending on N and A, which, as well as Uq and R, are considered finite. In the limits 
R ^ 0, Uq ^ oo, according to above given estimations, we have A ^ oo. 

Thus, when transforming to the singular potential the bound state energy remains finite (and, moreover, 
any given) value only if N will tend to infinity as 



3/2 _ 12V2an ( 



as R^Q. 
Let us write Ceup as 



K(0) 



Ceup = Ce ^ , ^1 2 , u ^^^"^ 



insert PH)) in ©and take account of the formulas 



dpiK(0)T4.(0) = irSn,k. (12) 



Then one obtains equations 



A ^ 



8^2^nu; £ / "^P'n + b + cpr ^^^^ 

n=0 



N 

N-2 \ ^ 



=0 
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Integrating the right side of Eq. (jl3p we obtain the foUowing transcendental equation 

N 

E 



8V27T ^\hujn-E\ 



(15) 



the roots of which determines the energy levels of electron in the considered combination of fields. This 
equation can be solved numerically (graphically); its roots a;„ — En/huj for different A is determined by 
the crossing of horizontal 8-\/2a7r/A with different values A with the function 



N 



1 



(16) 



It is seen from Eq. p6p that for attractive potential the main contributions in the sum gives small 
lengths Ax, which adjoins to the corresponding root with n to the left. Plot of f{x) as a function of x 
near the roots with n = 1,2, 3, 6, 10 is given in Fig. 1 for = 10®. 




FIG. 1: Plot of the function f{x) 

It follows from Eq. IjlSp that if the potential energy is small compared with huj then the energy level 
En is located between two energy levels hujn and huj{n + 1) with the exception of the negative energy 
level £^0- 

The energy levels can be approximately calculated if 



A 



87raV2 



< 1. 



Then, representing 



E„ 



TiLon + 5n, (5 < 



(17) 



for any number n = 0, 1, . . ., we obtain 



5n = 



327r2 



(18) 



First, let the singular potential can be considered as a perturbation therefore the coupling constant A < a. 
All the levels (5„ are shifted down with respect to the corresponding Landau levels hum for any n = 0, 1, . . .. 
In this case |^„| ^ Tluj. If £^0 < is the electron energy in the pure singular \5{y) potential then (5o is the 
magnetic field correction to the level Eq. 
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If the coupling constant A obeys the inequahty 8tt^/2 3> A > a then \E\ > Hlu and E is determined by 
Eq. p8| . It should be emphasized there exists the only energy level {Eq < 0) in a single attractive A^(r) 
potential. 

Plot of f{x) as a function of x near the root with n — 0, which corresponds ground (bound) state is 
given in Fig. 2 for N = 10®. 













-2035 












-2030 












-2025 












-2020 












-2015 












-2010 












-2005 
-2000 






-0.8 


-0.6 


X 


-0.4 


-0.2 






FIG. 2: Function f{x) near the root n = 

Graphically the root Xq = Eq/Tiuj is determined by the crossing of horizontal 8V2a7r/A with the hmction 
plotted in Fig. 2. We also see that the energy of bound electron state is analytical in A. 

The electron energy in bound state is not analytical in A for the model under discussion in the two- 
dimensional case. Indeed, Eg. p^ should be replaced by the equation 



N 



A 



47r n 

n=0 



6' 



E 

TujJ 



(19) 



for the two-dimensional case (see, Q)- Here we took account of all remarks concerning the number N, 
which were given when deriving of Eq. (jl5p . The energy of ground state (that is a negative root of 
Eg. p^ Eq < 0) can be easily found for ^ 1. Replacing the summation over n by the integration, we 
find as a result 



i.Ai^f^ 

47r \ b 



(20) 



and for the electron energy of bound state we obtain 



47r\ 

Eq = —hajNciq) ( j ■ 



(21) 



It is of importance to emphasize that in the two-dimensional case the coupling constant A is the 
dimensionless constant. Nevertheless, there exists the bound state in the attractive XS{r) potential. In 
the limits N oo, A — > we must require that N should depend on the dimensionless constant A so as 
the binding energy —Eq would remain finite as — > oo. Thus, the cutoff dimensionless parameter TV, 
which tends to infinity, transmutes in arbitrary binding energy |£^o|- This is the nonrelativistic analog 
the dimensional transmutation phenomenon. For the model under discussion in the absence of magnetic 
field this phenomenon was first considered in Ref. ^T§\ . 
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III. ELECTRON WAVE FUNCTION AND PROBABILITY CURRENT IN THE BOUND 

STATE 



Now we shall construct the wave function of electron in the negative-energy state with n — 0. It can 
be found from Eq.([71) by putting n — 0. Then, using Eqs. fTUl) . HH), we integrate over pi by means of the 
following integral^"' 



dxe Un{x + z)Uk{x + u) = - exp 

a 



iy{z + u) M + y 

h i[n — k) arctan ■ 



z — u 



Ink{p), (22) 



where the function Ink{p) = cxp(— p/2)L'j (p) is the Laguerre function of the argument 

+ {u- zf 

P= 2 ' ^ ' 

■^n~^{p) is the Laguerre polynomial and in our problem n = k ~ 0. 
Integration over p3 



'pI + 2m\Eo\ 



1= I dp3 . , (24) 



can be carried out in the complex plane, closing the integration contour into the lower half-plane for z > 
and in the upper half-plane for z < 0. The case A ^ a is of physical interest. For this case, one obtains 



^y/2m\E„\\z\/h 

, (25) 
V2rn\Eo\) 

Then, simple calculations leads to the normalized electron wave function in the form (see, also Refs. 

*o(r) - -^=- exp ^—^ ^ exp - V ' ^ ^ (26) 



2nla V J \ I J ^/2m\EQ 

It is of great interest that there is the probability current even when the electron is in the bound 
state in a superposition of a constant uniform magnetic field and a single attractive XS{r) potential. The 
probability current density is 

n ( x^ + y^ 2V2\z\\ ^ ( x' + y^ 2V2\z\\ 
= -2^^^''''' [-^^ — ) " 1^"^^^ — ) ' 

/ x^ + y^ 2V2\z\ \ ^ ^ 
Jy = jQxexp{ — 1 I' J^ = ^- (27) 

We see from Eq. (P7|) that the divergence of the probability current density is equal to zero everywhere. 
Therefore, the probability density 1^*0 (r)p and the current density ((27)) satisfy the continuity equation 

at ox oy oz at 

everywhere and the function |^'o(r)p is conserved in time. 

The electron wave functions in an external electromagnetic field are known to have the ambiguity, which 
is related to the ambiguity of the 4-potential (<&, A) of electromagnetic field. The latter is determined 
just with the exactness up to the gauge transformation 

A^A + VA, (29) 

c at 
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where / is an arbitrary function of coordinates. The Schrodinger equation does not change provided the 
replacement of the vector potential in the Hamiltonian is carried out simultaneously with the replacement 
of wave function accordingly 



1G 

*exp(-/ 



(30) 



In the considered case, we can cancel the phase factor in wave function ([26|) by means of the gauge 
transformation of vector potential with function / — Bxy/2. Under such a transformation the vector 
potential transforms from Eq. ([2]) to 



A== y(-2/,a;,0), 

the wave function (|26p becomes real but the (electric) current density 

being by the gauge invariant quantity, does not change. Note that 

J = K^*-^ 



(31) 



(32) 



(33) 



at the gauge ([?T|) . 

The probability current density is shown in Fig. 3 in which the coordinates x, y, z are measured in 
units of a. 




FIG. 3: Probability current density in space 



The vector field J can be treated semiclassically. Let the two-dimensional vector field j be a complex 
quantity 3 = jx + ijy in any plane z — constant whose the components jx and jy are functions of the 
complex variable v — x + iy and depend on z. The vector field 3{x,y) is a vortex field. Computing 
C = [V X j] , we obtain 



[V X j] = 



dx 



djx 
dy 



A 

TT 



VV 

2^ 



exp 



2a2 



(34) 



where A = ehF (z) / Ima^ , the function 



F{z) — exp — 



2\z\ 
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should be calculated on the plane z = constant and v* is the complex conjugate of v. In each plane 
z = constant the vector j{x,y) determines the vector field of a point-like vortex located at the point 
x,y = 0. 

One can write j in the form 



il{v) 
2Trv* 



(35) 



where 



liv) — Aw* exp —7-^ 



(36) 



is the vortex intensity that is the vector field circulation j (x, y) on any closed contour encircling the 
vortex at the point x,y = 0. 

The normalized two-dimensional vector field j(a;, y) is shown in the plane 2: = in Fig. 4 in which the 
coordinates x,y are measured in units of a. 
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FIG. 4: Probability current density in the 2 = plane 



One can see from Eq. (P7)l that the probability current density resembles a stack of "pancake" vortices" 
whose circulating "currents" around the magnetic field direction (z-axes) are mostly concentrated within 
the plane z ~ Q. 



IV. EFFECT OF HOMOGENEOUS ELECTRIC FIELD ON THE BOUND ELECTRON 

Now we consider the effect of a weak constant homogeneous electric field on an electron bound by an 
attractive singular potential in the presence of a constant uniform magnetic field. Related problems are 
the ionization of negative charged ions in accelerators as well as the removal of electrons from trapping 
potentials by a constant electric field. It is worthwhile to note that the problem concerning the removal of 
a charged nonrelativistic particle from a sp herically symmetric potential well by a constant homogeneous 
electric field was first partly solved in [l3| and the more rigorous formulas was given in ^J^] for the total 
and in [l^ for the differential probability. Here we consider this problem for the model studied with the 
inclusion of a weak constant homogeneous electric field specified by the potential 

U{z) = -\e\ez, (37) 

where e > is the electric field strength. 
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The inclusion of the electric field in the model under discussion leads to the quasistationary of the 
bound electron state and the appearance of the nonzero particle flux at large distances from the potential 
well. It is apparent that for considered electromagnetic field combination we need calculate only the 
particle flux in the z direction far from the potential well. 

Let us flnd the decay probability of the bound state per unit time, which is equal to the electron flux 
across the plane perpendicular to the z axis. To flnd this flux we need the Schrodinger equation solutions 
for an electron far from the potential well. So far as the electron flux is of interest in the z direction the 
magnetic held cannot be taken into account far from the well. Solutions of the Schrodinger equation for 
an electron in fleld ([57)1 are expressed via the Airy functions (see, for example, [l9j). 

So far as the particle flux has to be nonzero in the z direction we must choose the solution with the 
complex Airy function for the real z: 

i^p,^j„^E = e-'^*'^e'^^'''+^^y^'^V{Z), (38) 

where 

Pi,P2 are the eigenvalues of the operator of generalized momentum and V{Z) is the complex Airy function 
for the real Z. For Z < 1 



2z3/2\ i f 2z^/2 



exp I ) + - exp 



(40) 



Note that only the keeping of exponentially small term in solution (j40[) will give rise to the nonzero 
flux. 

Outside the well the electron wave function for large z can be written as the superposition of solutions 
in electric field with E = —Eq 

^(r) = e-'-E«*/'"' / dpidp2e'^P'''+P''y'^/''ViZ)F{pip2) (41) 



The appearance probability of electron far from the well in unit time w is equal to the flux across the 
plane perpendicular to the z axes 

The function F{pip2) is found for large z from Eq. (|4ip by the Fourier transform 

F{pip2) = ^ J^^.^. I dxdyV(r)e*(^°*-fi^~f^^)/'*. (43) 



47r2y(Z) 

In weak electric field e <C r7i^/2|_EoP^^/?i|e| = eo the electron wave function (|41]) differs insignificantly 



from the wave function of bound state ([26ll (we just consider this case) for \z\ ^ y ?i^eo/2m|£'o|e, so for 

such z we can substitute the function (|26ll instead of ^{r) in Eq. (|43l). 
Carrying out simple calculations, we obtain 

where C — l/(a-\/27r/) is the normalization constant of the wave function of bound state. Finally, for the 
total probability one obtains 



2y/T:Tia'^\e\e f ml 



(45) 
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V. RESUME 

It is of interest to compare wave function (j26p with the electron wave function of bound state in the 
only singular attractive h^\5{v)/2m potential. The latter can be easily obtained in the form 



where Iq = y ?i^/2m|£'o| and E'o < is the electron energy and the wave function (H^ is normalized as 
follows 

j \^{v)\^dV^l. (47) 

The distribution of probabilities of different coordinates of the electron 

I T , M , 9 , „ / 2r\ dr 
|*(r)|47rr'^dr = 2 exp 



/ '0 



is spherically symmetrical. Here r — -^/x^ + j/^ + z^. One sees that without magnetic field the electron 
is localized in the region ^ Iq and the probability current density for the state (|46)) is equal to zero 
everywhere. 

The distribution of probabilities of different coordinates of the electron in the state (^51 

a \ 2a^ L I a L 



has the cylindrical symmetry. Here p = ^/x"^ + y"^, z are the cylindrical coordinates. The electron is 
located in the region x,y ^ y/2a in the x, y plane perpendicular to the magnetic field and z ^ 21 along 
the magnetic field direction. 

It is of importance that the probability current arises even when the electron is in the bound state a 
superposition of a constant uniform magnetic field and a single attractive X5{r) potential. This (electric) 
current is the gauge invariant (physically observable) quantity and, therefore, it must involve (and exercise 
influence on) the electron scattering. 
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